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FOREWORD 

This document presents a final report of a research pro- 
gram (Contract NAS8-21485) performed by Lockheed's Huntsville 
Research & Engineering Center (Lockheed/Huntsville), for NASA- 
Marshall Space Flight Center (NASA/MSFC). 

Technical supervisor and coordinator of the contract were 
R. S. Ryan and H. J. Buchanan, respectively, both of the Aero- 
Astrodynamics Laboratory, NASA/MSFC. 
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SUMMARY 

A new approach which employs the combined system modes instead of 
the bending and slosh modes separately in formulating the mathematical 
model of a liquid propellant vehicle is used to simulate the atmospheric 
flights of the Saturn V vehicle. The model includes two rigid-body motions, 
four system bending modes and a control system with high-order control 
dynamics. The aerodynamic forces, thrust, gravitational field and vehicle 
vibrational characteristics which are functions of flight time are represented 
by 2 1 time- varying coefficients. Propellant slosh forces and slosh ampli- 
tudes may be calculated after the responses of the vehicle are obtained. 

An associated hybrid computer program was developed. The program 
uses the analog computer to perform the numerical integration of the system 
and the digital computer to compute and store the response parameters of 
the system when subjected to a large ensemble of wind samples. At the end 
of the first passage of the program, the mean values of the response param- 
eters are calculated. The statistics of the wind and vehicle response are 
computed in the second passage of the program. 

For an ensemble of 970 tape winds, the program takes less than five 
minutes of EAI 8900 hybrid computer time to provide the above data for 
Saturn V during its first 90 seconds of flight. The excellent computational 
speed suggests that the Monte Carlo hybrid simulation method is an effective 
and economic technique for conducting preliminary vehicle design studies. 


iii 
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Symbols 

a 

o 



c 


4i 


M^(25, t) 
M^(90, t) 
M B (25,t) 

Mg (90, t) 

Mp (25, t) 


NOMENCLATURE 


attitude feedback control coefficient, rad/rad 

attitude rate feedback control coefficient, rad/ rad/sec 

attitude coefficient associated with angle of attack, 
l/ sec^ 

attitude coefficient associated with engine deflection 
angle, l/sec^ 

til 

attitude coefficient associated with the i n system 

2 

bending mode, l/sec 

i** 1 system bending coefficient associated with engine 
deflection angle, m/sec^ 
til 

i system bending coefficient associated with angle of 
attack, m/sec^ 

lateral motion coefficient associated with attitude, 

/ 2 
m/ sec 

lateral motion coefficient associated with angle of 
attack, m/sec^ 

lateral motion coefficient associated with engine deflec- 
tion angle, m/sec^ 

bending moment coefficient due to aerodynamic force 
at station 2 5, N-m/ rad 

bending moment coefficient due to aerodynamic force 
at station 90, N-m/rad 

bending moment at station 2 5, N-m 
bending moment at station 90, N-m 

bending moment coefficient due to engine deflection at 
station 25, N-m/rad 


v 
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Symbols 
M' ' (90, t) 

|M\.(25,t)|. 




m 


ki 


t . 
ni 


V 


m 


V 


w 


< Y G>i 


< Y R>i 


NOMENCLATURE (Continued) 


bending moment coefficient due to engine deflection at 
station 90, N-m/rad 

tb 

i system bending moment coefficient at station 25, 
N-m/ m/ sec^ 

til 

i system bending moment coefficient at station 90, 
N-m/ m/ sec^ 

■fVi 

rotational acceleration coefficient associated with n 
th 

beam and i system bending mode, dimensionless 

slosh coefficient associated with k** 1 tank and i** 1 system 
bending mode, dimensionless 

til 

lateral acceleration coefficient associated with n beam 
th 

and i system bending mode, dimensionless 
vehicle velocity, m/sec 

wind velocity, m/sec 

th 

slope at gyro point due to the i system bending mode, 

l/m 

th 

slope at rate gyro due to the i system bending mode, 
l/m 

vehicle lateral displacement, m 


Greek 


(a 2>„ 

P E 


angle of attack, rad 

n"* 1 beam lateral acceleration, m/sec 
engine deflection angle, rad 


vi 
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Greek 




e 

n 



00 . 

1 


NOMENCLATURE (Continued) 


th 

slosh amplitude at k tank wall, m 

til 

damping ratio of the i system bending mode, 
dimensionless 

i** 1 system bending mode coordinate, m 
n t ^ 1 beam rotational acceleration, rad/sec^ 

til 

first slosh mode coordinate associated with k tank, m 
i** 1 system bending mode frequency, rad/ sec 


vii 
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Section 1 
INTRODUCTION 

To simulate the atmospheric flights of a liquid propellant space vehicle, 
the dynamics of the propellant are traditionally represented by a mass-spring- 
dashpot system (Refs. 1, 2 and 3). In the case of Saturn V, the mathematical 
model of the vehicle is limited to a dynamic system which consists of about 
10 second-order ordinary differential equations. The improsed restriction 
is due to the capacity of a hybrid computer or the excessive computer time 
which may be required for a digital simulation. Consequently, only the dy- 
namics of the propellant in three of the six tanks of the vehicle can be taken 
into account. The rest of the propellant is assumed to be solid masses 
attached to the vehicle. The model represents the vehicle adequately during 
liftoff; however, in the second half of the first- stage flight, propellant levels 
will decrease in height, and a proper description of the dynamic behavior of 
the propellant may require that higher slosh modes be considered. Further- 
more, the propellant which was taken as solid masses may have a larger 
share in the contribution to the dynamics of the vehicle than the sloshing 
propellant. 

Lockheed/Huntsville has introduced a system mode approach in form- 
ulating the mathematical model of a large complex liquid propellant space 
vehicle. The method uses the Lagrange equation for solving a coupled elas- 
tic and fluid system which characterizes a liquid propellant vehicle (Ref. 4). 
The modeling of a vehicle now considers the system bending modes instead 
of the bending modes and slosh modes of a vehicle. As a consequence, the 
dynamics of the liquid propellant will always be taken into account and the 
size of the mathematical model will be smaller than that obtained from the 
conventional approach. 
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After the mathematical model of a vehicle is adopted, the statistical 
response of the vehicle can be obtained analytically as well as numerically 
(Refs. 2, 3 and 5). Also, depending on the required accuracy, computation 
time, capability of the computer and the complexity of the model, the atmo- 
spheric flights of a vehicle can be simulated on an analog, digital or hybrid 
computer. Presently, the following four methods may be used for flight 
simulations (Ref. 6): 

1. Monte Carlo hybrid simulation, 

2. Monte Carlo digital simulation, 

3. State Space method and 

4. Sampling Model method. 


Although a judicious choice among the methods is controversial, the speed 
of the analog computation in performing the repetitive numerical integrations 
of a dynamic system suggests that hybrid simulation is definitely a good and 
economical method for conducting the preliminary design studies of a vehicle. 
The objective of this research effort is to show the feasibility of using the 
system mode approach in formulating mathematical models and of using 
hybrid flight simulations for Saturn V or future non-beamlike space vehicles. 


2 
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Section 2 

MONTE CARLO SIMULATION METHOD 

To ensure a successful launch of a space vehicle, the vehicle must be 
designed for all possible wind disturbances encountered during its atmo- 
spheric flight. Since the wind field is of a random nature and is difficult to 
define, a statistical approach should be used to determine the dynamic re- 
sponse of a vehicle to wind disturbances. One of the statistical approaches 

which is currently used in flight simulations is called the Monte Carlo 
simulation method (Ref. 6). The method employs the fast computational 
speed and unlimited storage spaces of modern computers to simulate the 
passage of a vehicle through a large number of recorded wind profiles. The 
statistics of the vehicle response are computed from the ensemble of the 
calculated responses. Depending on the required accuracy, computational 
time and the complexity of a mathematical model, the method can be executed 
on a digital computer or a hybrid computer. 

In general* the digital simulation provides good accuracy and, virtually 
no limitation on the size of a mathematical model. However, its application 
to dynamic statistical response studies is often prohibited by the excessive 
computer time which is needed in performing the repetitive numerical integra- 
tions of a system. In case an additional system mode or a nonlinear mathe- 
matical model needs to be considered, computational time will increase 
drastically. Fortunately, in preliminary design or parametric studies, the 
accuracy of the vehicle response becomes less critical if the simulation can 
show the basic characteristics of the system properly. In addition, the com- 
putational speed of an analog computer will not be affected if nonlinear mathe- 
matical models are considered. Hence, hybrid simulation may be used to 
obtain the statistical response of a dynamic system. Since the size of a hybrid 
simulation program often needs to be reduced to suit the capacity of a hybrid 


3 
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compute r, engineering judgment will play a major role in the modeling of a 
complex vehicle. The conclusion of this study is that digital simulation may 
be employed for spot checks of a dynamics analysis. Hybrid simulation should 
be used to conduct preliminary vehicle design studies. 


4 
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Section 3 

MATHEMATICAL MODEL 

The system mode approach is used to model the Saturn V vehicle. An 
associated hybrid computer program was developed to compute the statisti- 
cal response of the vehicle to 970 tape wind disturbances. In order to pro- 
vide a good quantitative study and, at the same time, to accommodate the 
hybrid computer program in the Lockheed/Huntsville computer, engineering 
judgment and compromises among accuracy, computational speed and other 
considerations were made in formulating the following mathematical model 
of Saturn V. The model includes two rigid-body motions, four system 
bending modes and a control system with high-order control dynamics. Data 
used in the simulation were obtained from different sources. Approxima- 
tions and extrapolations were also employed to define some of the vehicle 
parameters. Data related to the aerodynamic forces, thrust, nominal tra- 
jectory, gravitational field and the wind tape were obtained from Ref. 3 and 
R. S. Ryan of the Dynamics & Control Division, Aero-Astrodynamics Labo- 
ratory, MSFC. Vehicle bending and sloshing data were generated from the 
Tandem Beam-Slosh Computer Program developed under this contract 
(Ref. 4). 

The equations of motion of the model are given below. Further dis- 
cussions of these equations can be found in Refs. 2 and 7. 

a = <i> - 

m m 

y = K x 0 + K 2 a+ K 3 (3 e (2) 


5 
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= Cj a+ C 2 p E + c 4i T| i 


(3) 


i = 1 


f,-. = -2 £. W . -Hi -CO. T| i + D.a4 c 5 .(3 e (4) 


(i = 1,2, 3,4) 


Notations used in Eqs. (1) through (4) are defined below. Numerical data of 
the coefficients are given in Figs. A- 1 through A- 9 and Table A- 1 of Appen- 
dix A. Wind samples which have been converted as functions of flight time 
are provided for the simulation through an analog wind tape. 


a = angle of attack 

3 e = engine deflection angle 

th 

= damping ratio of the i n system bending mode 
th. 

r|^ = i system bending mode coordinate 

th 

= i system bending mode frequency 
<j> = attitude angle 

= attitude coefficient associated with angle of attack 

C 2 = attitude coefficient associated with engine deflection angle 

th 

c 4i = attitude coefficient associated with the i system bending 
mode 

th 

c 5i = * s Y stern bending coefficient associated with engine 
deflection angle 

th 

= i system bending coefficient associated with angle of attack 
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- lateral motion coefficient associated with attitude 

= lateral motion coefficient associated with angle of attack 

= lateral motion coefficient associated with engine deflection 
J angle 

V m = vehicle velocity 

V w = wind velocity 

y = vehicle lateral displacement 

The control system of the model is 


a o F($) 


-,) . T1 • 
r i 


4 

E <*': 

i = 1 


where 
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where s denotes the Laplace operator and 

p = 0.60051 
*0 

= 0.055535 
p 2 = 3.7209 x 10"' 

p* = 0.99253 
p* = 0.014998 
p* = 0.011931 
p* = 0.099515 x 10 
p^ = 0.021996 x 10 

Functions T(|3 ), F(4>) and F(<£) are the engine command, attitude and attitude 
rate transfer functions of the system, respectively. The coefficients a^, 

(Yy i and (Y» ) are the attitude gain, attitude rate gain, slope at gyro point 
associated with i^ system bending mode and slope at rate gyro associated 
with i** 1 system bending mode, respectively. Their numerical values are 
given in Table A- 1 of Appendix A. The above control system was modified 
by Ryan. Detailed discussion and additional references on the control system 
can be found in Ref, 8. 

The bending moments at stations 25 and 90 of Saturn V are calculated 
from Eqs. (6) and (7), respectively (Ref, 2). 


q =1.0 

q x = 0.24335 

q 2 = 8.2079 x 10“ 3 

q 3 = 0.090539 x 10' 
* 

q = 1.0 

q* = 0.62648 

q* = 0.12300 

q* = 8.7553 x 10" 3 

q* = 0.20157 x 10~ 3 

q* = 1.9072 x 10" 6 
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M (2 5, t) = M' (25, t) a + M' (25, t) p 
-d a p E 


+ 



(6) 


M n (90,t) = M' (90, t) a + M* (90, t) p 
J3 a p -Ci 


4 




Ml- (90, t) 



(7) 


whe re 


M 1 (25, t), M' (90, t) 
a a 


bending moment coefficients due to aerodynamic 
force at stations 25 and 90, respectively. 


M' p (25,t), M' (90,t) = 


bending moment coefficients due to engine deflec- 
tion at stations 25 and 90, respectively. 


ML (2 5, t) . , 

L *1 


ML (90, t) . 
r] ' 7 1 


th. 

i system bending moment coefficients 
at stations 25 and 90, respectively. 


Their values are presented in Figs. A- 10, A- 11 and Table A-l of Appendix A. 


The free vibrational characteristics of Saturn V are found by modeling 
the vehicle as a series of four beams (Fig. 1). The deflection of each beam 
is described by four fundamental deflection functions. The dynamics of the 
liquid propellant contained in each tank is represented by its first slosh mode. 
Bending properties of the vehicle are obtained from the NASA/MSFC Stodola 
tape. Slosh modes and frequencies of the liquid propellant contained in the 
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Saturn V tanks are calculated from the Lomen program (Ref. 9). Conse- 
quently, the generalized coordinates of this system are five beam end dis- 
placements, 16 fundamental beam deflection coefficients and six slosh 
coefficients. 

The system modes, frequencies and other quantities of the above 
Saturn Y model are computed from a Lockheed/Huntsville tandem beam 
slosh program, developed under the same contract. 

The first four system bending modes of Saturn V during its first 90 
seconds of flight are selected as those shown in Figs. A- 15 through A- 20 of 
Appendix A. 

Slosh amplitudes and local accelerations are calculated from Eqs. (8) 
through (11) after the dynamic response of the vehicle is found. 

4 

= b k Z) s ki ”i 

i = 1 

fa 2 ) n = -y+ Y. ‘ni'ii 

i = 1 
4 

0 =-</j + y p . ri. 

n ni 'i 

- i = l 

4 

^kl~ 2Z 4 S ki 

i = 1 


(8) 

(9) 

( 10 ) 
( 11 ) 
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where 


bj = 0.5229, 
b 4 = 0,6042, 


b 2 = 0.5390, 
b c = 0.6020, 

D 


b 3 = 0.5431, 
b 6 = 0.6630 


Notations used in Eqs, (8) through (11) are defined below. Numerical values 
of the coefficients are given in Figs. A- 12 to A- 14 of Appendix A. 


6^ = slosh amplitude at tank wall 

(a 0 ) = beam lateral acceleration 

v 2 n 

= beam rotational acceleration 


n 


th 


kl 


ki 


t . 
m 


ni 


= first slosh mode coordinate associated with k tank 

= slosh coefficient associated with k^ tank and i ^ system 
bending mode 

= lateral acceleration coefficients associated with n*'* 1 
beam and i^ system bending mode 

= rotational acceleration coefficient associated with n 
beam and i^ system bending mode 


th 


With the above information, the lateral slosh force distribution along the 
vehicle longitudinal axis can be calculated. Detailed discussions and deriv- 
ations of the system mode approach are given in Ref. 4. 


12 
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Section 4 

HYBRID COMPUTER PROGRAM 

Lockheed/Huntsville has an inhouse EAI 8900 hybrid computer which 
provides digital, analog and hybrid capabilities for solving complex engi- 
neering and scientific problems. The digital computer is general purpose 
with a high-speed 16K memory and a complete complement of software. 

The analog computer performs parallel integration to obtain rapid solutions 
to systems of differential equations in real time or faster. Special software 
is used to program the hybrid interface, set up and check out simulation, 
and debug programs. 

A hybrid program based on the mathematical model defined in Section 
3 was developed to simulate the atmospheric flights of Saturn V. The pro- 
gram uses the analog computer to perform the repetitive numerical integra- 
tion of the system and the digital computer to control program logic while 
calculating and storing the response parameters of the system. In order to 
accommodate the model with the currently available equipment of the com- 
puter, curve fitting was employed to approximate some of the time -varying 
coefficients. The program is designed to simulate the first 90 seconds of 
the flight and to record the response parameters of the vehicle at every five 
seconds of flight time without slowing down the analog computation. At the 
end of the first passage of the program, the mean values of the vehicle re- 
sponses are computed. Standard deviations of the response parameters are 
computed in the second passage of the program. The number of exceedance 
counts, envelopes and plots of certain response parameters are obtained 
during the run. Consequently, the program will provide the following features 
in a single run: 

1. Response of an individual flight 

2. Envelope of a response parameter 
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3 . Mean values and variances of response parameters 

4. Number of exceedance counts for a response parameter 

5. Statistics of the wind. 

For an ensemble of 970 tape winds, the program takes less than five 
minutes of EAI 8900 hybrid time to compute the above data for Saturn V 
during its first 90 seconds of flight. The excellent execution speed obtained 
in this study justifies the conclusion that the Monte Carlo hybrid simulation 
technique is indeed an effective method for determining the statistical re- 
sponse of a launch vehicle. 

The wiring diagrams of the analog program are given in Appendix B. 

A listing of the digital program is presented in Appendix C. Detailed dis- 
cussions and additional references on hybrid computation can be found in 
Refs. 10 and 1 1. 


14 
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Section 5 

EXAMPLE - STATISTICAL RESPONSE OF SATURN Y 

The mathematical model defined in Section 3 was used to obtain the 
response of the Saturn V vehicle to wind disturbances* Effort was made to 
show that the system mode approach and hybrid simulation are good methods 
for determining the statistical response of a liquid propellant space vehicle. 
Additional parametric studies and improvement of the model are needed in 
order to provide a good quantitative response study of Saturn Y. 

The developed hybrid program provides the following information: 

1. Traces of fig, a, y, M B (90,t), M B (25,t), V w> rjj, p 2 » r) 3 > t]^, 

(j> , F (<$>), (j) , F (<t>) and timing pulses; 

2. Number of exceedance counts of p^, M B (90,t) and M B (25, t) through 
comparators; 

3. Number of exceedance counts of 6^, 6^ and 6^ through digital com- 
putation; 

4. The mean values, and 3a of R-,, (90, t), (25, t) and V , 

sLt r5 ±5 W 

where a denotes the standard deviation of a variable; 

5. The mean values of a, <j> , 8, , (a-,) , 6 and £,, , where k = 1,2, ... , 

ic Lt n n ic i 

6 and n = 1,2,3. 

As an example, a model which includes the first three system bending modes 
was used to compute the statistics of Saturn V to 970 tape winds. Output of 
the hybrid program for the first 90 sec of flight was studied. Some of the 
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results are presented below. The envelope of the 970 tape winds is shown in 
Fig. 2. The NASA/MSFC synthetic wind profile (Ref. 3) was used to check 
out the program. Traces of the variables during this slow run are given in 
Fig. 3. Fig. 4 provides " stretched- out" traces of the variables during fast 
runs to tape wind disturbances. The scales of these traces are defined in 
Table A- 2 of Appendix A. The mean values and 3cr values of V^, |3g, Mg 
(90, t) and Mg(25,t) are given in Fig. 5. Percentage of the number of ex- 
ceedance counts of Mg(90,t), Mg (25,5), 6^. (k = 4, 5, 6) are given in Fig. 6. 


16 
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Fig. 2 - Envelope of Tape Winds 
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w 
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Fig. 3 - Traces of Variables of the Checkout Run 
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Fig, 4 - Stretched Traces of the Variables During Fast Runs 
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Fig. 4 - (Continued) 
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Section 6 

CONCLUSIONS AND RECOMMENDATIONS 

The advantage of using the system mode approach in formulating a 
mathematical model of Saturn V for its atmospheric flight simulations was 
demonstrated. Future space vehicles, such as space shuttles, will have non- 
beamlike configurations. Based on the present preliminary designs, the 
vehicle will carry a large amount of liquid propellant during its atmospheric 
flight. In addition, a large number of structural modes must be included to 
model the vehicle properly. The system mode approach is possibly the only 
method which could provide a good mathematical model for flight simulations 
on a computer. Also, this research effort has shown that the Monte Carlo 
hybrid simulation technique is an excellent method in obtaining the statistical 
response of a vehicle. It is an especially powerful method for conducting 
preliminary design or parametric studies of a vehicle. 

With the experience which has been obtained from various flight sim- 
ulation studies of Saturn V (Ref. 6), one could carry the present state of art 
one step further. For instance, additional research effort may be made in 
the following areas. 

1. The Sampling Rate of Hybrid Simulations : The wind tapes which 

Lockheed/Huntsville and Computer Sciences Corporation (Ref. 3) used in the 
hybrid and analog simulations, respectively, are functions of flight time. 
This means that the atmospheric flights of a vehicle are simulated on a pre- 
determined trajectory. It would be preferable for the wind tape to be a 
function of altitude. Then, flight simulations could be perturbed on any 
arbitrary nominal trajectory with the necessity of making only one wind tape. 
Furthermore, if a nonlinear or an advanced mathematical model should be 
used, the trajectory of a flight may be defined as a parameter which depends 
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on both time and vehicle response. However, two problem areas need to be 
resolved before the new wind tape can be used in a hybrid simulation: 

• The wind data may have to be stored in the digital computer 
and fed into the analog computer at discrete points during a 
computer run. Delays due to the repetitive ADC and DAC 
operations need to be studied. 

® The sampling rate of a hybrid simulation must be defined 
such that the high frequency components of the wind can be 
properly presented. 

2. Hybrid Simulation of Nonlinear Mathematical Models: Techniques 

and programs developed for flight simulations have all been based on linear 
mathematical models. The modeling of a vehicle often leads to non- 
linear dynamic systems. Since the hybrid simulation is so fast, and non- 
linear systems present no difficulty in analog computation, accurate solutions 
may be obtained if nonlinear mathematical models can be also used in flight 
simulations. 

3. Sampling Model Method for Flight Simulations: The sampling 

model method is a compromise between the methods of State Space and 
Monte Carlo simulation (Ref. 5). In this method, a sampling model picks 
winds from the wind sample population in accordance with its appropriate 
probability distribution. The same statistics of the vehicle response may 
be obtained by simulating the passage of a vehicle through only part of the 
wind samples. Since each wind sample has its own time history, to define 

a sampling model may be more complicated than solving the problems which 
can be treated directly by applying the standard techniques of statistics 
(Refs. 12 and 13). However, this approach may eliminate the limitations 
which are often encountered in other simulation techniques, such as accu- 
racy, capacity and execution time. 
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Appendix A 

VEHICLE PARAMETERS AND COEFFICIENTS 
OF SATURN V 
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Fig. A- 12 - Slosh Amplitude Coefficient 
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Fig. A-12- (Continued) 
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Fig. A- 12 - (Continued) 
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Fig. A- 12 - (Continued) 
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Fig. A- 14 - (Continued) 
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Fig. A- 15 - System Bending Modes at t=0 
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Fig. A- 15 - (Continued) 
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Fig, A- 16 - System Bending Modes at t=20 sec 


A-22 


LMSC/HREC D149410 


SA'Soa ^-qsam moocl 


5 t.nt 

6£N. MASS =3.®I/S40* 



A-23 




LMSC/HREC D149410 


&A~M)3 MOOSt T z 20.0 

WOOS (is ) 

Fscaur&scv s 2.329 

♦ 03 

6£W. MASS 33.S5SX1Q 



Fig. A- 16 - (Continued) 
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Fig. A- 17 - System Bending Modes at t=40 sec 
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Fig. A- 17 - (Continued) 
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Fig* A- 17 - (Continued) 
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Fig* A- 18 - System Bending Modes at t=60 sec 
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Fig. A- 18 - (Continued) 


A-32 



LMSC/HREC D1 49410 


SA- JO 3 4-0g A M 1*00 ZL 
MOD 202) 

a d./S$ 

&SM. WAS 3 *S.S4®3I0^ C 




ill 





Ml 


8 CL A T I VC AMPLITUDE 


4th System 
Bending Mode 
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Fig. A- 19 - System Bending Modes at t = 80 sec 
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Fig. A-20 - System Bending Modes at t = 100 sec 
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Table A- 2 

VARIABLE SCALES 


Variable (s) 

Unit/ Division 

3e ’ 

0.00125 

a 

0.025 

© • 
y 

0.12 

M b (90, t) 

4 x 10 4 

M B (25,t) 

10 6 

V 

2 

w 


’ll 

0.02 

’l 2 

0.01 

’l 3 

0.005 

y ) 4 

0.003 

<t> , F ($) 

0.002 

i . f (4>) 

0.001 
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Appendix B 

ANALOG WIRING DIAGRAMS 













Fig. B- 1 
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Fig. B-3 - Analog Wiring Diagram (Attitude Filters) 
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Fig. B-5 - Analog Wiring Diagram (Wind, Angle of Attack and Altitude) 
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Fig. B-7 - Analog Wiring Diagram (Coefficients) 
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Fig. B-10 - Analog Wiring Diagram (Logic for Wind Tape and A/D) 
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LISTINGS OF THE DIGITAL PROGRAM 
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SJOB 

SFR. . IU1 » *4003 

LOGICAL INDIC 

extended title 

REAL M90X « M25X * M90B « M25B* M90 V ♦ M25V « M90 * M25 * MA25 * M A90 * MB90 « M903 • 
IM253 ♦ I 18* I2B* I3B 
DAT A I OCT 8X0 00377/ 

COMMON / BLOCK / I OCT 8*1*1 ERR * NM AX * NRM AX « NP ASS ♦ NRUN * INDIC* I CHAN ♦ 

1 Tl * T2 *N » I LOC 

COMMON ADSCALt 16> *DASCALC 16) * BX ( 4 ) * DEL4X ( 3 ) « DEL5X ( 3 ) «DEL6X(3> * 

1 IBX(4) ♦ I M90X ( 6 ) ♦ IM25X <6) * 104(3) « 105(3) ♦ ID6( 3 ) *BB( 18 ) * AB ( 18 ) * 

2YDDB( 18) *PHI2B( 18 > *VWB< 18) *ATA1B( 18) « ATA2B ( 18) «ATA3B( 18) «ATA4B( 18) 
3* 

4ATD 1 B ( 18) » ATD2B ( 18) • ATD3B< 18) * ATD4B< 1 8) *PH IB< 18 ) « XI 1B< 18 ) *X 1 28 ( 18 ) 
5«XI3B( 18 ) «XI4B( 18) *Xl5B( 18) »XI6B( 18) « DEL 1 B < 18) «DEL2B( 18) «DEL3B( 18) 
6»DEL4B( 18) *DEL5B( 18) *DEL6B( 18) « BV<18)« 

7AV (18)* VWV ( 1 8 > * DAC ( 8 > , ADC (16)*B(18)«A(18)» YDD ( 1 8 ) « PH I 2 ( 1 8 ) * PH I ( 1 8 ) 
8* AT A 1 ( 18) » ATA2< 18) ♦ ATa3< 18) « ATA4< 18) » ATD 1 (18) *ATD2( 18) «ATD3( 18) ♦ 
9ATD4 ( 18) *VW( 18) *DEL4( 18) *DEL5( 18) *DEL6( 1 8 ) « I I BX ( 4 ) » I IM90X(6) * 

1 I IM25X ( 6 ) ♦ I ID4(3) * I ID5(3)« I ID6(3) *B3( 18) * A3C 18) *VW3( 18) * TITLE (20) ♦ 
2S11 ( 18) * S 1 2 ( 18) »S13( 18) »S14( 18) « S2 1(18) »S22( 18) * S23 ( 18) «S24( 18) « 
3S31 (18) *S32( 18) *S33( 18) «S34( 18) »S41 ( 18) *S42( IS) «S43( 18) ,S44( 18) « 
4S51 (18) *552(18) *S53( 18) »S54( 18) *S61 ( 18) «S62( 18) «S63( 18) * S64 ( 18) ♦ 
5Tli(18)«Tl2(18)»T13(18)»Tl4(18)*T21( 18 ) ♦ T22 ( 1 8 ) * T23 ( 18) *T24< 18) * 
6T31 (18) *T32( 18) *T33( 18) *T34( 18)*Pll(l8)*P12(18)«Pl3(18)«Pl4(18). 
7P21 (18) *P22( 18) *P23( 18> »P24( 18) »P31 ( 18) *P32< 18) *P33( 18) «P34( 18) 

1 ♦ A21B( 18) * A22B( 18) *4236(18) 

COMMON 

1 M90X(6) « M25X ( 6 ) *M90B( 18) »M25B( 18) *M90V( 18) «M25V( 18) «M90( 18) * 

1 M25 ( 1 8 ) ♦ MA25 ( 1 8 > ♦ MA90 ( 1 8 ) « MB90 ( 1 8 ) * M903 ( 1 8 ) ♦ M253 ( 1 8 ) « 

21 1B( 18) ♦ I2B( 18) * I3B( 18) 

DAT A 1 0CT8/0 00377/ VRAN0040 

1 000 FORMAT ( 2 1H»C0NS0LE* DISCONNECTED) 

1010 FORMAT! 29H* NON-EX I STING COMPONENT REQ • D ) 

1030 FORMAT ( 8 I 1 0 ) 

1040 FORMAT ( 1 8H OVERFLOW DAC CH. 12) 

1050 FORMAT (27H NON-EXISTING DAC CH. REQ • D ) 

1060 FORMAT ( 24H MISSED BEGINNING OF RUN) 
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1070 FORMAT ! 26H ERROR IN BLOCK ADDRESSING) 

1080 FORMAT! 18H OVERLOAD ADC CH. 12) 

1090 FORMAT! 27H NON-ExISTInG ADC CH* REQ»D) 

1100 FORMAT! 14H END OF PASS 1) 

1110 FORMAT! 20A4) 

1150 FORMAT! 1X9E14.7/) 

1160 FORMAT! 1X91 14 ) 

1700 FORMAT! I 14.8E14.7) 

C INPUT DATA 

c input a/d and d/a scale factors 

14 READ 1 1 10* (TITLE! 1 > ♦ 1=1 *20) 

CALL READ 

C INITIALIZE ANALOG MODE AND START OF RUN INDICATOR 
lOO CALL Q I CO! l * I ERR) 

IF ! IERR*EQ* 1 > GO TO 5 
IF! IERR*EQ*2 ) TYPE 1000 
PAUSE 
GO TO 1 00 
5 CONTINUE 

C INPUT NO. OF POINTS/RUN AnD NO. OF RUNS 
READ 10?0*NMAX*NRMAX 
C SET INITIAL CONDITIONS 
21 NPASS=1 

I BX ( 1 )sO 

I BX ! 3 ) = 0 
DO 8 1=1 *5*2 
IM90X! I ) =0 
IM25X! I ) =0 
8 CONTINUE 

DO 230 1=1*3 
I D4 ! I ) = 0 
ID5! I )=0 
I 06! I >=0 
230 CONTINUE 

DO 9 1=1.18 
BB ! I ) =0. 

AB! I ) =0 * 

VWB! 1 )=0, 
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ataibc I ) = o. 

ATA2B ( I ) *0. 

ATA3B( I ) *0. 

AT A4B ( I ) =0* 

ATD 1 B ( I ) =0. 

ATD2BI I ) =0® 

ATD3B=0 • 

ATD4B® 0 • 

YDDB( I )=0. 

PH I 2B ( I 5=0. 

PHIBC I )*0. 

9 CONTINUE 

C INITIAL CONDITION FOR PASS 2 

I BX < 2 ) =0 
IBX< 4 ) *0 

40 DO 16 I =2*6*2 
IM90XC I > =0* 

I M25X ( I > =0 * 

16 CONTINUE 

DO IT 1=1*18 
BV( I ) =0 « 

AV< I >=0. 

M90V ( I ) = 0. 

M25V ( 1 ) = 0. 

VWV< I 5=0. 

17 CONTINUE 
4 NRUN=1 

CALL QTEFF0( 1 *0* INDIC* IERR) 

IF( IERR*EQ. 1 ) GO TO 51 
I F ( I ERR*EQ*2 ) TYPE 1000 
I F ( I ERR • EQ • 3 ) TYPE 1010 
PAUSE 
GO TO 4 

C TEST FOR PASS NO** 

51 I F < NPASS.EQ. 1 ) GO TO 10 
I F ( NPASS.EQ* 2 ) GO TO 1 1 
C SET FIRST SET OF EXCEEDANCE LEVELS 
10 DAC( 1 )=-BX( 1 )*DASCAL( 1 ) 
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DAC ( 2 ) =-BX ( 3 ) *DASCAL ( 2 ) 

DAC C 3 ) =-M9QX C 1 ) *DASCA|_ < 3 ) 

DAC ( 4 ) =-M90X { 3 ) *DASCA(_ < 4 > 

DAC ( 5 ) =~M90X ( 5 > *DASCAl * 5 ) 

DAC ( 6) = -M25X C 1 ) *DASCA|_ ( 6 ) 

DAC « 7 ) =“M25X ( 3 ) *DASCAL ( 7 ) 

DAC ( 8 ) * -M25X ( 5 ) *DASC AL < S ) 

GO TO 12 

C SET SECOND SET OF EXCEEDANCE LEVELS 

11 DAC( 1 )=~BX(2)»DASCAL( l > 

DAC ( 2 ) =-BX < 4 ) *DASCAL ( 2 > 

DAC < 3 ) JS-M90X < 2 ) *DASCAl < 3 ) 

DAC ( 4 ) »-M90X < 4 ) *DASCAl < 4 ) 

DAC ( 5 ) *-M90x ( 6 ) *DASCAl < 5 ) 

DAC ( 6 ) =— M25X ( 2 ) *DASCAl< 6 ) 

DAC ( 7 ) =-M25X ( 4 ) *DASCAl ( 7 ) 

DAC ( 8 ) =-M25X < 6 ) *DASCAl < 8 ) 

12 CALL QDAJM0(0*8«DAC« lERR* 1CHAN) 
I F < 1 ERR • EQ • 1 ) GO TO 3 

1F< IERR*EQ.2)TYPE 1040* I CHAN 
IF< 1ERR»EQ»3 ) TYPE 1050 
PAUSE 
GO TO 14 

C TEST FOR BEGINNING OF RUN 

3 CALL QTEFFO ( 1 «0* IND1C ♦ I ERR) 

IF< IERR.EQ.1 ) GO TO 13 
IF( IERR*EQ.2)TYPE 1000 
IF( IERR.EQ.3) TYPE IOjO 
PAUSE 
GO TO 14 

13 IF< INDIC) GO TO 18 

C TEST FOR BEGINNING OF RUN 

19 CALL QTEFF0( 1 «0« INDIC* IERR) 

I F ( INDICJGO TO 20 
GO TO 19 

20 CALL QTEFFO (1*0*1 ND I C ♦ I ERR 1 
IF< INDICJGO TO 20 
GO TO 50 
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C MISSED BEGINNING OF RUN 
18 TYPE 1060 

CALL QIC0( 1 « I ERR) 

NRUN=NRUN-1 
PRINT i 1 60 *NRUN 
PRINT 1 1 50 iTl 
PRINT 1150, T2 
PAUSE 

C TEST FOR PASS NO* 

IF(NPASS*EQ*1 )GO TO 21 
IF<NPASS*EQ*2 )GO TO 40 
C SELECT OPERATE TO START ANALOG 
50 CONTINUE 
/ LDOB =*2047**12 

CALL CLOCKS 
IF< IERR*EQ* 1 )GO TO 22 
IF( IERR*EQ*2 ) TYPE lOOo 
TYPE 1060 
PAUSE 

C READ ADC'S 

22 DO 23 N s 1 « NMAX 

CALL QSTREO ( 0 * 16* I ERR ) 

I F ( I ERR • EO • 1 ) GO TO 24 
IF( IERR.EQ.2)TYPE 1070 
CALL Q I CO ( 1 * IERR) 

PAUSE 
GO TO 14 

24 CALL QADCVOtO* 16*ADC* IERR* ICHAN) 
I F ( I ERR • EQ • 1 ) GO TO 41 
IF( IERR *EQ*2 ) TYPE 1080 * ICHAN 
IF( IERR*EQ*3 ) TYPE 1090 
CALL Q I CO ( 1 » IERR) 

PAUSE 
GO TO 14 

41 B<N)=ADC< 1 )*ADSCAL< 1 ) 

A <N) = ADC C 2 ) *ADSCAL ( 2 ) 
YDD(N)=ADC(3)*ADSCAL( 3> 

PH I 2 ( N ) * ADC ( 4 ) *ADSCAL ( 4 ) 
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PHI ( N ) * ADC ( 5 ) *ADSCAL ( 5 > 

ATA1 ( N ) = ADC ( 6 ) *ADSCAL ( 6 ) 

ATA2 { N ) * ADC ( 7 ) *ADSCAL ( 7 > 

ATA3(N)« ADC ( 8 ) * ADSC AL ( 8 ) 
ATA4(N)=ADC<9)*ADSCAL(9) 

ATD1 ( N) *ADC( 10)*ADSCAl( 10) 

ATD2 ( N) * ADC < 1 1 >*ADSCAl< 1 1 > 

ATD3(N)*ADC( 12)*ADSCAl< 12) 

ATD4 ( N ) * ADC (13) *ADSCAl ( 13) 

VW( N ) * ADC 1 14)*ADSCAL( 14) 

M90 ( N ) * ADC (15) *ADSCAL (15) 

M25 ( N > * ADC (16) *ADSC AL (16) 

23 CONTINUE 

C SELECT HOLD MODE TO STOP ANALOG 
/ LDOB * * 2447 * « 12 

C READ EXCEEDANCE COUNT 

CALL QLGINO ( l ♦ ILOC* I ERR ) 

I F ( IERR*EQ* 1 ) GO TO 25 
IF< IERR»EQ*2 ) TYPE lOOO 
PAUSE 
GO TO 14 

C SELECT THE IC MODE ON THE ANALOG 
25 CONTINUE 

/ LDOB * * 3047 • ♦ 1 2 

T1*CL0CKR(0) 

C TEST PASS NO. 

IF ( NPASS.EQ. 1 ) GO TO 26 
IF ( NPASS.EQ*2 ) GO TO 2 

C CALCULATE THE FIRST SET OF EXCEEDANCE COUNTS 
C** 

C* SUBROUTINE COUNTS 

C DO 27 1*1 *5.2 

C 26 I IBX( 1 ) =0 

C I IBX( 3) *0 

C I I M90X ( I )=0 

C I I M25X ( I ) *0 

C 27 CONTINUE 

C I F ( ILOC .AND. 1)11 BX ( 1 ) = 1 
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C I F ( I LOC • AND • 2 ) I I BX ( 3 ) « 1 

C !F( IL0C*AND*4) I IM90X< i )=1 

C IF ( I LOC* AND *8 > I I M90X ( 3 ) = 1 

C I F ( I LOC • AND* 16)11 M90X ( 5 ) = 1 

C I F I I LOC* AND *32 ) I I M25X ( 1 ) a 1 

C IF ( ! LOC» AND»64 ) I I M25X (3 ) * 1 

C I F ( I LOC • AND • 1 28 ) ! I M25X ( 5 ) = 1 

C I BX ( 1 ) * I BX { 1 ) + 1 I BX < 1 ) 

C I BX ( 3 ) » I BX ( 3 ) + 1 I BX ( 3 ) 

C DO 28 I=l*5«2 

C I M90X ( I ) = I M90X ( I ) + 1 1 M90X ( I ) 

C I M25X ( I ) = I M25X ( I ) + IIM25XU) 

C 28 CONTINUE 
C DO 210 1=1 .3 

C I ID4( I ) *0 

C I ID5< I ) *0 

C I 1D6( 1) =0 

C 210 continue 

C DO 200 1=1 «NMAX 

C DEL4( I )=.6042*<S41 ( I )*ATA1 ( I )+S42< I )*ATA2( I )+S43( I )*ATA3( I > + 

c 1 S44C I )*ATA4C I ) ) 

C DEL5( I >=«59*(S51 ( I )*ATA1 ( I >+S52< I >*ATA2< I )+S53( I )*ATA3< I ) 

C 1 + S54 { I )*ATA4< I ) ) 

C OEL6 ( I )**64*<S6i ( I )*ATA1 < I >+S62< I )*ATA2< I )+S63( I )*ATA3( I )+ 

c 1 S64 < I )*AtA4< I ) ) 

C I F C DEL4 ( I ) «GT*DEL4X< 1 ) ) I ID4( 1 > = 1 

C I F < DEL4 < I > *GT*DEL4X<2) >1 ID4( 2) « 1 

C I F ( DEL4 ( I ) *GT«DEL4X(3) ) I ID4(3)=1 

C I F ( DEL5 ( I ) «GT*DEL5X( 1 ) ) I ID5( 1 )=1 

C I F ( DEL5 < I ) •GT®DEl_5X<2} ) I ID5( 2 ) * 1 

C IF ( DEL5 ( I ) • GT * DEL5X < 3 > ) I ID5< 3) = 1 

C IF( DEL6( I ) #GT *DEL6X< 1 ) ) I ID6( 1 ) = 1 

C IF(DEL6C I ) •GT*DEI_6XC2) ) I ID6( 2) = 1 

C IF ( DEL6( I ) *GT *DEL6X(3) ) I ID6(3)=1 

C 200 CONTINUE 
C DO 220 1=1 *3 

C I D4 ( I ) = I D4 ( I ) + I I D4 ( I ) 

C I D5 C I ) = I D5 < I ) + I I 05 < I ) 
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C ID6( I ) » ID6( I )+I ID6( I > 

C 220 CONTINUE 

C CALCULATE THE SUMMATION OF XI 
C DO 29 1*1 *NMAX 

C BB( I ) »BBC I )+ B ( I ) 

C AB( I >=AB< I )+ At I ) 

C M90B( I >=M90Bt I ) +M90 ( I ) 

C M25B< I )=M25B( I 1+M25C I > 

C VWB ( I ) * VWB t I ) +VW ( I ) 

C ATAlBt I )*ATA1B< I ) +ATA 1 < I ) 

C ATA2B< I )=ATA2B( I )+ATA2< 1 ) 

C ATA3B< I >=ATA3BC I )+ ATa3<I> 

C ATA4B( I >*ATA4B< I )+ ATo4 t I) 

C ATDIBt I )*ATD1B( I J+ATDl ( I ) 

C ATD2B( I >*ATD2Bt I )+ATD2t I ) 

C ATD3B< I >=ATD3B< I )+ATD3( I ) 

C ATD4B( I >=ATD4B( I )+ATD4< I ) 

C PHIB< 1 )*PHIB( I )+ PHI ( I > 

C YDDB ( I ) = YDDB < I ) +YDD ( I ) 

C PHI2B< I >*PHI2B< I )+PHl2< I ) 

C 29 CONTINUE 
C RETURN 

C END 

C** 

26 CALL COUNTS 

IFt NRUN*EQ«NRMAX >GO To 30 
T2*CLOCKR(0 ) 

NRUN=NRUN+1 
GO TO 3 

C CALCULATE MEAN VALUES 
30 DO 31 I * 1 ♦ NMAX 

BBt I > =BB( I ) /NRMAX 
AB< I ) =AB( I l/NRMAX 
M90B« I )=M90B< I ) /NRMAX 
M25B ( I )=M25B( I ) /NRMAX 
PHIBt I ) * PH I B t I ) /NRMAX 
VWB ( I ) »VWB ( I ) /NRMAX 
ATAlBt I )= ATAlBt 1 ) /NRMAX 
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ATA2B ( I ) = ATA2B ( I ) /NRMAX 
ATA3B< I )«ATA3B( I ) /NRMAX 
ATA4B< I )*ATA4B( I >/NRMaX 
ATDIBt ! )=ATD1B< I > /NRMAX 
ATD2B( 1 >=ATD2B( 1 > /NRMAX 
ATD3B( I )*ATD3B( I ) /NRMAX 
AT04B ( I )*ATD4B( I J/NRMAX 
VDDB ( I > * YDDB < I ) /NRMAX 
PH 1 2B ( I )=rRHI2B< I J/NRMAX 

31 CONTINUE 
C END OF PASS 1 

TYPE 1 1 00 
PAUSE 

NPASS=NPASS+1 
GO TO 4 

C CALCULATE THE SECOND SET OF EXCEEDANCE COUNTS 
2 I |BX<2)»0 

I I BX ( 4 ) =0 
DO 32 I=2.6«2 
I IM90XC I > =0 
I IM25XC I ) *0 

32 CONTINUE 

IF< I LOC • AND* l ) I IBX<2)=1 
IF< ILOC*AND*2> I IBX(4) =1 
IF < ILOC.AND*4) I IM90X(2) = 1 
IF < I LOC* AND * 8 ) I I M90X( 4 ) = 1 
IF ( I LOC • AND • 1 6 ) I IM90X'6>*1 
IF< ILOC* AND»32) 1 IM25X(2)*l 
IF ( ILOC* AND.64 ) I IM25X(4)=1 
IF ( I LOC* AND. 128) I IM25X<6)=1 
ASSEMBL 


CA 

/I LOC 

AHM 1 

/I IBX+l .*15 

ROT 

1 

AHM 1 

/IIBX+3..15 

ROT 

1 

AHM 1 

/I IM90X+1 * . 15 

ROT 

1 
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AHM1 

/ I IM90X+3* • 15 

ROT 

1 

AHM1 

/l IM90X+5* ♦ 15 

ROT 

1 

AHM1 

/l IM25X+1 Mlg 

ROT 

1 

AHM1 

/ I IM25X+3* » is 

ROT 

1 

AHM1 

/ 1 I M25X+5 * • 1 5 


fortran 

I BX ( 2 ) = I BX < 2 ) + 1 1 BX ( 2 ) 

I BX ( 4 > s 1 BX C 4 ) + 1 1 BX ( 4 ) 

DO 33 I =2* 6*2 

IM90XC 1 ) * I M90X ( 1 ) + I 1 M90X ( 1 ) 

IM25XC 1 >*IM25X< I )+I IM2SXC I ) 

33 CONTINUE 

C CALCULATE SUMMAT I ON ( X-XB ) **2 
DO 34 I = I » NMAX 

BV< I )*BV( I )+ ( B ( I ) — BB ( I ) ) **Z 
M90 V ( I ) = M90V < I ) + ( M90 ( I > -M90B ( I ) ) **2 
M25V ( I ) =M25V ( I) + < M25 ( I ) -M25B ( I ) ) **2 
VWV< I )*VWV( I ) + ( VW ( I ) -VWB ( I ) ) **2 

34 CONTINUE 

IF(NRUN*EQ.NRMAX)GO To 35 
NRUN=NRUN+1 
GO TO 3 
C* 

C SUBROUTINE VAR I AN 

C CALCULATE THE VARIANCES 
C 35 DO 42 I * I ♦ NMAX 
C BV< I )=BV( I ) / ( NRMAX— 1 ) 

C A V ( I ) a AV ( I ) / ( NRMAX- 1 ) 

C M90 V ( I ) =M90V ( I )/( NRMAX- 1 ) 

C M25V< I )=M25V< I) / (NRMAX- 1 ) 

C VWV < I > aVWV ( I )/( NRMAX- i ) 

C B3C I )=3**SQRT(BV( I ) ) 

C A3 ( I )=3.*SQRT( AVC I ) ) 

C M903( I )=3**SQRT(M90V( j ) ) 
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C M253 < I ) = 3 • *SQRT ( M25V ( | ) ) 

C VW3< I ) »3 • *SQRT ( VWV ( I ) ) 

C XI 1B( I >s SI 1 ( I >*ATA1B( I ) +S12( I )*ATA2B< I ) +S13< I )*ATA3B< I ) 

C 1 +SI4C I»*ATA4B( I ) 

C X 1 2B ( I )* S21 ( I )*ATA1B( I ) +S22( I )*ATA2B( I ) +S23( I >*ATA3B< I ) 

C 1 +S24 < I )*ATA4B ( I ) 

C X I3B ( I ) * S3I < I )*ATA1B( I ) +S32 ( I ) *ATA2B ( I ) +S33 ( I ) *ATA3B ( I) 

C I +S34 ( I >*ATA4B( I ) 

C XI4B< I )= S41 ( I ) *ATA 1 B ( I ) +S42 < I) *ATA2B ( I ) +S43 < I ) *ATA3B < I) 

C 1 +S44( I )*ATA4B< I ) 

C X I5BC I ) * S51 ( I ) *AT A 1 B ( I ) +S52 ( I ) *ATA2B< I ) +S53 ( I ) *ATA3B < I > 

C 1 +S54< I )*ATA4B< I ) 

C I 1 B ( I ) *— PHI2B( I ) +Pl l < I )*ATD1B< I )+P12< I )*ATD2B< I )+Pl3< I )*ATD3B< I ) 

C I +P14< I)*ATD4B( I ) 

C 1 + T34< I )*ATD4B< I ) 

C 1 +P24< I >*ATD4B< I ) 

C I 2B ( I ) =-PHI2B< I ) +P21 < I )*ATD1B( I )+P22< I )*ATD2B< I )+P23< I )*ATD3B< I ) 

C X I6B( I ) = S61 ( I )*ATA1B< I) +S62< I )*ATA2B( I ) +S63( I ) *ATA3B( I ) 

C 1 +S64< I )*ATA4B( I ) 

C DEL 1 B( I ) = XI 1B( I >*.5299 

C DEL2B ( I )= XI2B<I)*,539 

C DEL3B( I >= XI3B< I )**543l 

C DEL4B < I > = XI4B( I )*«6042 

C DEL5B< I )= X I SB ( I )*.59 

C DEL6BI I ) = X I 6B ( I )**64 

C A21B( I ) = YDDB < I ) + T 1 1 ( I ) *ATD 1 B ( I > +T \ 2 < I ) *ATD2B ( I ) +T13( I )*ATD3B< I ) 

C 1 + T 14 ( I )*ATD4B< I ) 

C A22B ( I ) * YDDB ( I > + T21 ( I ) *ATD 1 B ( I ) +T22 ( I ) *ATD2B ( I) +T23< I )*ATD3B< I ) 

C 1 + T14( I )*ATD4B( I ) 

C A23B ( I ) * YDDB ( I ) + T3 1 ( I ) *ATD 1 B ( I ) +T32 ( I ) *ATD2B ( I ) +T33< I )*ATD3B< I ) 

C 1 3B ( I ) = —PH I 2B < I ) +P3 1 ( I >*ATDIB< I )+P32( I )*ATD2B( I > +P33 C I )*ATD3BC I ) 

C 1 +P34 ( I )*ATD4B< I ) 

C 42 CONTINUE 
C RETURN 

C END 

C* 

35 CALL VAR I AN 
CALL WRITE 
PAUSE 
GO TO 14 
END 
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SUBROUTINE WRITE 

REAL M90X » M25X * M90B » M25B*M90V« M25V < M90 « M25 » MA25 *MA90 ♦ MB90 * M903 » 

1 M253 ♦ I 1 B« I 2B ♦ I3B 

COMMON ADSCAL { 16) *DASCAL( 16) ♦ BX ( 4 ) *DEL4X(3) *DEL5X(3) «DEL6X(3> ♦ 

1 IBX(4) % I M90X ( 6 ) * I M25X ( 6 ) ♦ ID4(3) ♦ 105(3) t 1 06 ( 3 ) « BB ( 18 ) ♦ AB ( 1 8 ) « 

2YDDB ( 18) »PHI2B( 18) »VWb( 18) ♦ ATA1B( 18) « ATA2B ( 18)* AT A3B ( 18) « ATA4B ( 18) 
3# 

4 ATD 1 B ( 1 8 ) * ATD2B (18)* ATD3B (18)* ATD4B (18)*PHIB(18)«XI1B(18)«XI2B(18) 
5»XI3B( 18) »XI4B( 18) *XI5B( 18) «XI6B( 18) * DEL 1 B ( 18) «DEL2B( 18 ) « 0EL3B ( 18) 
6«DEL4B( 18) *DEL5B( 18) »DEL6B ( 18) * BV ( 18)* 

7AV( 18) *VWV( 18) *DAC(8) *ADC( 16)«B(18)*A( 18) *YDD( 18) *PHI2( 18) «PHI < 18 ) 
8 » AT A 1 ( 18) *ATA2( 18) *ATa3( 18) *ATA4( 18) « ATD 1 ( 18) «ATD2( 18) «ATD3( 18) * 
9ATD4( 18) *VW( 18) *DEL4( 18) *DEL5( 18) *DEL6( 18) • I IBX(4) * I IM90X(6) « 

1 I IM25X( 6) * I 104(3) ♦ I I D5( 3 ) * I I D6( 3 ) *B3 ( 18 ) * A3 ( 1 8 ) *VW3( 18) *TITLE(20) ♦ 
2S11(18)*S12(18)*S13(18)*S14(18)«S21(18) «S 22 < 18) *S23( 18) «S24( 18) . 
3S31 (18) * S32 (18)* S33 ( 18) *S34( 18) *S41 ( 18) *S42( 18) tS43( 18) «S44( 18 ) « 
4S51 (18) *S52( 18) *S53( 18) *S54( 18) *S61 ( 18) * S62 ( 1 8 ) « S63 ( 18) *S64( 18 ) « 
5Tll(18)*Tl2(18)*T13(l8)*T14(18)tT21(18) *T22( 18) «T23( 18) *T24( 18 ) . 
6T31 (18) * T32 (18) ♦ T33 ( 18) *T34< 18)»P1K18)«P12(18)*P13(18)*P14(18)« 
7P21 ( 18) *P22( 18) *P23( 18) *P24( 18) *P3l ( 18) *P32( 18) *P33( 18) »P34( 18) 

1 »A21B( 18) »A22B( 18) ♦ A23B( 18) 

COMMON 

1 M90X(6) *M25X(6) *M90B( 18) *M25B( 18) *M90V( 18) *M25V( 18) «M90( 18) « 

1M25 ( 18) « MA25 ( 18) *MA90< 18) *MB90( 18) *M903( 18) «M253( 18) « 

21 1B( 18) * I2B( 18) * I3B( 18) 

1120 FORMAT ( 1 HI *5X*20A4 ) 

1130 FORMAT ( 1 HO* BOH THE EXCEEDANCE LEVELS ARE LISTED ABOVE THE COUNTS) 

1 140 FORMAT ( 1 HO * 5HBX < 1 ) 9X5HBX ( 2 ) 9X5HBX ( 3 ) 9X5HBX ( 4 ) ) 

1150 FORMAT ( 1 X9E1 4 *7/ ) 

1 160 FORMAT ( 1X9114) 

1 170 FORMAT ( 1H0»8HMB90X< 1 ) 6X8HMB90X ( 2) 6X8HMB90X ( 3 ) 6X8HMB90X ( 4 ) 6X 
1 8HMB90X ( 5 ) 6X8HMB90X ( 6 ) ) 

1 180 FORMAT ( 1 HO » 8HM825X ( 1 > 6X8HMB25X ( 2 ) 6X8HMB25X ( 3 ) 6X8HMB25X ( 4 ) 6X 
1 8HMB25X ( 5 ) 6X8HMB25X ( 6 ) ) 

1190 FORMAT ( 1H0.71HTHE VARIABLES ARE LISTED ABOVE THE MEAN SIGMA SQUARE 
1 * AND 3 SIGMA VALUES) 

1200 FORMAT ( 1H0»4HB( 1 ) 10X4HB(2) 10X4HB(3) 1 0X4HB ( 4 ) 10X4HB(5) 10X4HB<6> 10X4 
1 HB ( 7 ) 1 0X4HB ( 8 ) * 1 0X4HB ( 9 ) ) 
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1210 FORMAT ( 1 HO ♦ 5HB ( 10)9X5hB( 1 1 )9X5HB( 12)9X5HB( 13)9X5HB( 14)9X5HB( 15)9X5 
1 HB ( 16) 9X5HB ( 17)9X5HB( 1 8 > » 

1220 FORMAT C 1H0,4HA( 1 ) |0X4HA(2) 10X4HA(3> 10X4HA(4) 10X4HA< 5) 1 0X4HA ( 6 ) 1 0X4 
1 HA <7)1 0X4HA (8)1 0X4HA ( g ) ) 

1 230 FORMAT (IHOi 5HA < 10) 9X5HA (11) 9X5HA (12) 9X5HA (13) 9X5HA ( 14) 9X5HA ( 15) 9X5 
1 HA <16) 9X5HA( 17)9X5HA( 18) ) 

1 240 FORMAT ( 1 HO » 7HMB90 ( 1 > 7X7HMB90 ( 2 ) 7X7HMB90 ( 3 ) 7X7HMB90 ( 4 ) 7X7HMB90 ( 5 ) 7X 
1 7HMB90 ( 6 ) 7X7HMB90 ( 7 ) 7X7HMB90 ( 8 ) 7X7HMB90 (9) ) 

1 250 FORMAT ( 1 HO * 8HMB90 (10) 6X8HMB90 (11) 6X8HMB90 (12) 6X8HMB90 (13) 6X8HMB90 ( 
114) 6X8HMB90 (15) 6X8HMB90 (16) 6X8HMB90 (17) 6X8HMB90 (18)) 

1260 FORMAT ( 1H0«7HMB25( 1 ) 7X7HMB25 ( 2 ) 7X7HMB25 ( 3 ) 7X7HMB25 ( 4 ) 7X7HMB25 ( 5 ) 7X 
1 7HMB25 ( 6 > 7X7HMB25 ( 7 ) 7X7HMB25 ( 8 ) 7X7HMB25 ( 9 ) ) 

1 270 FORMAT ( 1 HO « 8HMB25 (10) 6X8HMB25 (11) 6X8HMB25 (12) 6X8HMB25 (13) 6X8HMB25 ( 
114) 6X8HMB25 (15) 6X8HMB25 ( 16) 6X8HMB25 (17) 6X8HMB25 (18)) 

1 280 FORMAT ( 1 HO « 5HVW ( 1 ) 9X5HVW ( 2 ) 9X5HVW ( 3 ) 9X5HVW ( 4 ) 9X5HVW ( 5 ) 9X5HVW (6 )9X 
1 5HVW ( 7 ) 9X5HVW ( 8 ) 9X5HVW ( 9 ) ) 

1290 FORMAT ( 1 HO »6HVW ( 10)8X6HVW( 1 1 )8X6HVW( 12)8X6HVW( 13)8X6HVW( 14)8X 
1 6HVW( 1 5 ) 8X6HVW ( 16)8X6HVW( 17)8X6HVW( 18) ) 

1300 FORMAT ( 1H0«46HTHE VARIABLES ARE LISTED ABOVE THE MEAN VALUES) 

1310 FORMAT ( 1 H0«6HPHI ( 1 ) 8X6HPH I <2)8X6HPHI ( 3 ) 8X6HPH I (4)8X6HPHI (5)8X 
1 6HPHI (6)8X6HPHI (7)8X6HPHI (8)8X6HPHI (9) ) 

1320 FORMAT ( 1 HI ) 

1330 FORMAT ( 1H0*7HPHI ( 10)7X7HPHI ( 1 1 ) 7X7HPHI ( 1 2 ) 7X7HPH I ( 1317X7HPHI ( 1 4 ) 7X 
1 7HPHI ( 15) 7X7HPH I ( 16) 7X7HPHI ( 17)7X7HPHI ( 18) ) 

1340 FORMAT ( 1H0«8HDEL4X( 1 ) 6X8HDEL4X ( 2 ) 6X8HDEL4X ( 3 ) ) 

1 34 1 FORMAT ( 1 HO * 7HDEL 1 ( 1 ) 7X7HDEL 1 ( 2 > 7X7HOEL 1(3) 7X7HDEL 1(4) 7X7HDEL 1 ( 5 > 7X 
1 7HDEL 1 ( 6 ) 7X7H0EL 1 (7) 7X7HDEL l (8) 7X7HDEL 1(9)) 

1342 FORMAT ( 1H0.8HDEL1 ( 10)6X8HDELI ( 1 1 )6X8HDEL1 < 12)6X8HDEL1 ( 13) 6X 

1 8HDEL1 ( 14) 6X8HDEL1 ( 15)6X8HDEL1 ( 1616X8HDEL1 ( 17)6X 

1 8HDEL 1 ( 18) ) 

1343 FORMAT! 1H0«7HDEL2( 1 ) 7X7HDEL2 ( 2 > 7X7HDEL2 ( 3 ) 7X7HDEL2 ( 4 ) 7X7HDEL2 ( 5 ) 7X 

1 7HDEL2 ( 6 ) 7X7HDEL2 ( 7 ) 7X7HDEL2 ( 8 ) 7X7HDEL2( 9 ) ) 

1350 FORMAT ( 1 HO »8HDEL5X ( 1 ) 6X8HDEL5X ( 2 > 6X8HDEL5X ( 3 > ) 

1360 FORMAT ( 1H0«8HDEL6X( 1 ) 6X8HDEL6X ( 2 ) 6X8HDEL6X ( 3 ) ) 

1 370 FORMAT ( 1 HO . 8HDEL2 (10) 6X8HDEL2 (11) 6X8HDEL2 (12) 6X8HDEL2 ( 1 3 ) 6X 
1 8HDEL2 (14) 6X8HDEL2 (15) 6X8HDEL2 (16) 6X8HDEL2 ( 1 7 ) 6X 

1 8HDEL2 ( 18) ) 

1 380 FORMAT ( 1H0. 7HDEL3 ( 1 ) 7X7HDEL3 ( 2 ) 7X7HDEL3 ( 3 ) 7X7HDEL3 ( 4 ) 7X7HDEL3 ( 5 ) 7X 
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1 

1390 FORMAT ( 
1 
1 

1400 FORMAT ( 
1 

1410 FORMAT ( 
1 
1 

1420 FORMAT ( 
1 

1430 FORMAT ( 
1 
t 

1440 FORMAT ( 
1 

1450 FORMAT ( 
1 
1 

1460 FORMAT ( 
1 
1 

1470 FORMAT ( 
1 
1 

1480 FORMAT < 
1 
1 

1490 FORMAT ( 
1 
1 

1500 FORMAT ( 
1 
1 

1510 FORMAT ( 
1 
1 

1520 FORMAT ( 


7HDEL3 ( 6 ) 7X7HDEL3 < 8 ) 7X7HDEL3 < 9 ) ) 

1H0»8HDEL3( 10)6X8HDEL3( 1 1 )6X8HDEL3( 12)6X8HDEL3< 13>6X 
8HDEL3C 14>6X8HDEL3( 15>6X8HDEL3< 16)6X8HDEL3( 17>6X 
8HDEL3< 18) ) 

1 HO » 7HDEL4 < 1 > 7X7HDEL4 ( 2 > 7X7HDEL4 ( 3 ) 7X7H0EL4 ( 4 ) 7X7HDEL4 < 5 ) 7X 
7HDEL4 ( 6 > 7X7HDEL4 < 7 ) 7X7HDEL4 < 8 ) 7X7HDEL4 ( 9 ) ) 

1 HO * 8HDEL4 (10) 6X8HDELA (11) 6X8HDEL4 (12) 6X8HDEL4 ( 1 3 ) 6X 
8HDEL4 (14) 6X8HDEL4 (15) 6X8HDEL4 (16) 6X8HDEL4 ( 1 7 ) 6X 
8HDEL4 ( 18) ) 

1 HO « 7HDEL5 ( 1 > 7X7HDEL5 ( 2 > 7X7HDEL5 ( 3 ) 7X7HDEL5 ( 4 ) 7X7HDEL5 ( 5 ) 7X 
7HDEL 5(6) 7X7HDEL5 ( 7 ) 7X7HDEL5 ( 8 ) 7X7HDEL5 ( 9 ) ) 

1 HO * 8HDEL5 (10) 6X8HDEL5 (ID 6X8HDEL5 (12) 6X8HDEL5 ( 1 3 ) 6X 
8HDEL5( 14) 6X8H0EL5 ( 15)6X8HDEL5( 16)6X8H0EL5( 17) 6X 
8H0EL5( 18) ) 

1 HO ♦ 7MDEL6 ( 1 ) 7X7HDEL6 ( 2 > 7X7H0EL6 ( 3 ) 7X7HDEL6 ( 4 ) 7X 7HDEL6 ( 5 ) 7X 
7HDEL6 ( 6 > 7X7HDEL6 ( 7 ) 7X7HDEL6 ( 8 ) 7X7HDEL6 ( 9 ) ) 

1 HO *8HDEL6( 10)6X8H0EL6( 1 1 >6X8H0EL6( 12)6X8HDEL6( 1316X 
8HDEL6( 14)6X8HDEL6( 15>6X8HDEL6( 16)6X8HDEL6( 17)6X 
8HDEL6( 18) ) 

1 HO * 8HALPH 1(1) 6X8HALPH 1 ( 2 ) 6X8HALPH 1(3) 6X8HALPH 1 ( 4 ) 6X 
8HALPH1 ( 5 ) 6X8HALPH1 ( 6 ) 6X8HALPH 1 (7)6X8HALPH1 (8) 6X 
8HALPH1 (9) ) 

1 HO » 9HALPH 1 ( 10) 5X9HALPH 1 ( 1 1 )5X9HALPH1 (12) 5X9HALPH 1 ( 13 )5X 
9HALPH1 ( 14) 5X9HALPH 1 ( 15)5X9HALPH1 ( 16) 5X9HALPH 1 ( 17 )5X 
9HALPH1 ( 18) ) 

1 HO * 8HALPH2 ( 1 ) 6X8HALPH2 ( 2 ) 6X8HALPH2 ( 3 ) 6X8HALPH2 ( 4 ) 6X 
8HALPH2 ( 5 ) 6X8HALPH2 ( 6 ) 6X8HALPH2 ( 7 ) 6X8HALPH2 ( 8 ) 6X 
8HALPH2(9) ) 

1 HO * 9HALPH2 (10) 5X9HALPH2 (11) 5X9HALPH2 (12) 5X9HALPH2 < 1 3 ) 6X 
9HALPH2 ( 1 4 ) 5X9HALPH2 ( 15 ) 5X9HALPH2 ( 16)5X9HALPH2( 17 )6X 
9HALPH2 ( 18) ) 

1 HO ♦ 8HALPH3 ( 1 ) 6X8HALPH3 ( 2 ) 6X8HALPH3 ( 3 ) 6X8HALPH3 ( 4 ) 6X 
8HALPH3 ( 5 ) 6X8HALPH3 ( 6 ) 6X8HALPH3 ( 7 ) 6X8HALPH3 ( 8 ) 6X 
8HALPH3 ( 9 ) ) 

1 HO . 9HALPH3 (10) 5X9HALPH3 (11) 5X9HALPH3 (12) 5X9HALPH3 ( 1 3 ) 5X 
9HALPH3 ( 14 )5X9HALPH3( 15 ) 5X9HALPH3 ( 1 6 ) 5X9HALPH3 ( 17)5X 
9HALPH3 ( 18) ) 

1 HO * 7H I DD 1 ( 1 ) 7X7HIDD1 (2)7X7HIDD1 (3)7X7HIDD1 ( 4)7X7HIDD 1 ( 5 ) 7X 
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1 

1530 FORMAT ( 
1 
1 

1540 FORMAT ( 
1 

1550 FORMAT < 
1 
I 

1560 FORM AT < 
1 

1570 FORMAT ( 
1 
1 

1580 FORMAT C 
1 

1590 FORMAT C 
1 

1600 FORMAT ( 

1 

1610 FORMAT C 

1 

1620 FORMAT < 

1 

1630 FORMAT < 

I 

1640 FORMAT « 

1 

1650 FORMAT ( 

1 

1660 FORMAT < 

1 

1670 FORMAT C 
I 

1680 FORMAT C 
I 

1690 FORMAT ( 

1 


7HIDD1 (6)7X7HIDD1 (7J7X7HI0D1 (8)7X7HIDD1(9) ) 

1H0.8HIDD1 ( 1 0 ) 6X8H IDD 1(11) 6X8HIDD 1(12) 6X8H I DD 1(13) 6X 
8H I OD 1(14) 6X8H I DD 1 C J 5 ) 6X8H I DD 1(16) 6X8H I DD 1 ( 1 7 ) 6X 
8HIDD1 ( 18) ) 

1H0,7HIDD2<1 > 7X7HI DD2 ( 2 ) 7X7H I DD2 ( 3 ) 7X7HIDD2 < 4 ) 7X7H I DD2 ( 5 ) 
7X7HIDD2 ( 6 ) 7X7HIDD2 ( 7 ) 7X7H IDD2 ( 8 ) 7X7HID02 ( 9 ) ) 

1 HO « 8H 1 002 ( 10>6X8HIOD2( 11 >6X8HlDD2( 12 > 6X8HIDD2 ( 1 3 ) 6X 
8HIOD2 ( 14 )6X8H1DD2 ( 15 ) 6X8HIDD2 ( 16 )6X8HIDD2( I 7>6X 
8HIDD2( 18) ) 

1 HO * 7H I D03 C 1 ) 7X7H I DD3 ( 2 ) 7X7H I DD3 ( 3 ) 7X7H I DO 3 ( 4 ) 7X7H I DD3 ( 5 ) 7X 

7H I DD3 ( 6 ) 7X7H I DD3 ( 7 ) 7X7H I DD3 ( 8 ) 7X7H I DD3 ( 9 ) ) 

1 HO * 8M 1 DD 3 ( 1 0 ) 6X8H1DD3 ( 1 1 )6X8H1DD3( 12)6X8HIDD3( 13)6X 
8HIDD3( 14)6X8HIDD3( 15)6X8HIDD3( 16)6X8HIDD3( 17J6X 
8HIDD3( 18) ) 

1H0.6HXI1 U)8X6HX11(2)8X6HXI1<3)8X6HXI1C4)8X6HXI1(5)8X 

6HXI 1 <6)8X6HXI 1 (7>8X6HXI 1 <8)8X6HXI 1 (9) ) 

1 HO ♦ 7HX 1 1 ( 1 0 ) 7X7HX 1 1 < 1 1 17X7HXI 1 ( 1217X7HXI 1 ( 13)7X7HXI1 ( 14 )7X 
7HXI 1 ( 15 ) 7X7HX 1 1 ( 16>7X7HXI 1 < 17J7X7HXI 1(18)) 

1 HO * 6HX I2<1) 8X6HX 12(2) 8X6HX 12(3) 8X6HX 12(4) 8X6HX 12(5) 8X 
6HX 12(6) 8X6HX 12(7) 8X6 HX 12(8) 8X6HX 12(9)) 

1 HO » 7HX 1 2 ( 10) 7X7HX 12(11) 7X7HX 1 2 ( 12) 7X7HX I 2 ( 13) 7X7HX I 2 ( 14) 7X 
7HX 1 2 ( 15) 7X7HX 1 2 ( 16) 7X7HX 1 2 ( 17) 7X7HX I 2 ( 18) > 

1 HO « 6HX 13(1) 8X6HX 13(2) 8X6HX 13(3) 8X6HX 13(4) 8X6HX 1 3 ( 5 ) 8X 
6HX 13(6) 8X6HX 13(7) 8X6HX 13(8) 8X6HX 13(9) ) 

1 HO » 7HX I 3 ( 10) 7X7HX 13(11) 7X7HX 1 3 ( 12) 7X7HX I 3 ( 13) 7X7HX 13 ( 14 ) 7x 
7HX 1 3 ( 15) 7X7HX I 3 ( 16) 7X7HX I 3 ( 17) 7X7HX I 3 ( 18) ) 

1 H0*6HX 14(1) 8X6HXI4 ( 2 ) 8X6HX I 4 ( 3)8X6HX 14(4) 8X6HX I 4 ( 5 )8X 
6HX 14(6) 8X6HX 14(7) 8X6HX 14(8 ) 8X6HX 14(9) ) 

1 HO « 7HX 14(10) 7X7HX I 4 ( 1 1 ) 7X7HX 1 4 ( 1 2 ) 7X7HXI 4 ( 1 3 ) 7X7HX 1 4 ( 1 4 ) 7X 
7HX I 4 ( 15) 7 X 7HX I 4 ( 1 6 ) 7X7HX 1 4 ( 17) 7X7HX 1 4 ( 18) ) 

1 HO • 6HX I 5 ( 1 )8X6HXI5(2)8X6HXI5(3)8X6HXI5(4)8X6HXI5(5)8X 

6HX 15(6) 8X6HX 15(7) 8X6HX 15(8) 8X6HX 15(9) ) 

1 HO ♦ 7HX I 5 ( 10 ) 7X7HX 1 5 ( 11) 7X7HX 1 5 ( 1 2 ) 7X7HX 1 5 ( 1 3 ) 7X7HX 1 5 ( 1 4 ) 7X 
7HX 15(15) 7X7HX 15(16) 7X7HX 15(17) 7X7HX 15(18)) 

1 HO * 6HX 16(1 ) 8X6HX 16(2) 8X6HX 16(3) 8X6HX 16(4) 8X6HX I 6 ( 5 ) 8X 
6HX 16(6) 8X6HX 16(7) 8X6HX 16(8) 8X6HX 16(9) ) 

1 HO « 7HX I 6 ( 10) 7X7HX 1 6 ( 1 1 > 7X7HX 16 ( 1 2 ) 7X7HX 1 6 ( 1 3 ) 7X7HX I 6 ( 1 4 ) 7X 
7HX 16(15) 7X7HX 16(16) 7X7HX 16(17) 7X7HX 16(18)) 
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1710 FORMAT ( 1H1 ) 

C PRINT EXCEEDANCE COUNTS 

PRINT 1 1 20 ♦ (TITLE! I > « 1= 1 « 20 ) 
PRINT 1130 
PRINT 1 140 

PRINT 1 150* (BX( I ) , 1=1 ,4) 
PRINT 1 160. ( IBX( I > » 1 = 1 .4) 
PRINT 1170 

PRINT 1 150. (M90X( I ) . Isl .6) 
PRINT 1 160. ( IM90XI I ) • I=1.6> 
PRINT 1 180 

PRINT 1150. ( M25X ( I ) . 1=1 .6) 
PRINT 1 160. ( IM25X( I ) ♦ 1=1 .6) 
PRINT 1340 

PRINT 1 150. ( DEL4X ( I ) . 1=1 .3) 
PRINT 1 160. ( ID4( I ) . 1=1 .3) 
PRINT 1350 

PRINT 1 ISO. (DEL5X( I > . 1=1 .3) 
PRINT 1 160. ( IDS( I ) . 1=1 .3) 
PRINT 1360 

PRINT 1 150. (OEL6X( I ) « 1= 1 .3) 
PRINT 1 160, ( ID6( I ) . 1=1*3) 

C PRINT MEAN VALUES AND VARIANCES 


PRINT 

1 190 







PRINT 

1200 







PRINT 

1 iso. 

(BB( I ) i 

> I 

X 

1 

♦ 9 

) 

PRINT 

1 iso. 

( BV ( I ) i 

► I 

X 

1 

♦ 9 

) 

PRINT 

1150, 

( B3 ( I ) i 

i I 

s 

1 

49 

) 

PRINT 

1210 







PRINT 

1150. 

(BB( I ) < 

» I 

s 

1 

04 

18) 

PRINT 

1 iso. 

<BV( I > i 

> I 

X 

1 

04 

18) 

PRINT 

1 150, 

(B3U ) i 

> I 

s 

1 

0 • 

18) 

PRINT 

1240 







PRINT 

1 150, 

<M90B< I > 

4 

I 

= 1 

49 ) 

PRINT 

1 150, 

( M90V ( I ) 

4 

I 

sl 

49 ) 

PRINT 

1150, 

( M903 ( I ) 

4 

I 

si 

49 ) 

PRINT 

1710 







PRINT 

1250 
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PRINT 1 150* (M90BC I ) ♦ I=lO« 18) 
PRINT 1 150» (M90V( I) • Isl0» 18) 
PRINT 1150* (M903 ( I ) * 1=1 0* 18) 
PRINT 1320 
PRINT 1260 

PRINT 1 150* (M258 ( I) » 1=1 *9) 
PRINT 1 1 50 • < M25 V (I)*I=1*9) 
PRINT 1 150« < M253 ( I > * 1 = 1 *9) 
PRINT 1270 

PRINT 1 1 50 « ( M25B (I)*I=lO*18) 
PRINT 1 ISO* (M25VI I ) « IalO* 18) 
PRINT l ISO. (M253C I ) * I=lO. 18) 
PRINT 1280 

PRINT I ISO* (VWB< I ) ♦ 1 = 1 *9) 
PRINT 1 150* (VWVt I ) * 1 = 1 *9) 
PRINT 1 ISO* (VW3( 1 ) ♦ I = l *9) 
PRINT 1290 

PRINT 1 ISO* ( VWB< I > « I = l0* 18) 
PRINT 1 150* (VWV( I ) * I=iO* 18) 
PRINT 1 ISO* (VW3U ) « I = lO* 18) 
PRINT 1300 
PRINT 1220 

PRINT 1 ISO. (ABC I )* 1 = 1 *9) 

PRINT 1230 

PRINT 1150. (ABCI )*I=10*18> 
PRINT 1310 

PRINT 1 150* (PHIBC I >* 1*1 *9) 
PRINT 1330 

PRINT 1150* (PHIBC I >» 1*10,18) 
PRINT 1710 
PRINT 1341 

PRINT 1 150* (DEL1BC I ) * 1 = 1 *9) 
PRINT 1342 

PRINT 1 150. (0EL1BC I ) ♦ 1=10* 18) 
PRINT 1343 

PRINT 1 150* (DEL2BC I) * 1 = 1 *9) 
PRINT 1370 

PRINT 1 150* (0EL2BC I ) ♦ 1*10* 18) 


C- 18 



LMSC/HREC D1 49410 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


1380 

1 150* <DEL3B< I > « 1=1*9) 


1390 

1 150* (DEL3B< I ) « 1=10* 18) 
1400 

1150* <DEL4B< I ) » 1*1 *9) 
1410 

1 150. <DEL4B< I > * 1=10. 18) 
1420 

1 150. <DEL5B<I> ♦ 1=1 *9) 


1430 

1 150. (DEL5B< I > * 1 = 10* 18) 


1440 

1150* (DEL6B< I ) * 1 = 1 .9) 
1450 


1150* (DEL6B« I ) * 1 = 10* 18) 
1710 


1460 


1 150* (A21B< I ) * 1*1*9) 
1470 

1 1 50 *<A21B(I)*lslO«18) 
1480 

1150,(A22B<I)*I=1.9) 

1490 

1 150* (A22BI I > *1=10*18) 
1500 

1 150* (A23B( 1 > *1=1*9) 

1510 

1 150* (A23B( 1 ) *1=10.18) 

1520 

1 ISO, { 1 1B( I ) * 1 = 1 *9) 

1530 

1 150, ( I 1B( I ) * 1 = 1 0. 18) 
1540 

1 150. ( I2B( I ) * 1 = 1 .9) 

1550 

1 150. ( I2B( I ) . I = lO* 18) 
1560 
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PRINT 1 150* < I3B< I ) * 1 = 1 *9> 
PRINT 1570 

PRINT 1 150* ( I3B( I) * I = l0» 18) 
PRINT 1710 
PRINT 1580 

PRINT 1 ISO* (XI 1B< l ) * 1*1 « 9 ) 
PRINT 1590 

PRINT 1 150* (XI 1B( I) • IstlO* 18) 
PRINT 1600 

PRINT 1 1 50 ♦ < X I 2B ( I ) *1=1*9) 
PRINT 1610 

PRINT 1 150t (XI2B( I ) • IslO* 18) 
PRINT 1620 

PRINT 1 150<<XI3B< I ) ♦ 1*1 *9) 
PRINT 1630 

PRINT 1 1 50 « < X I 3B (I)*I=lO*18) 
PRINT 1640 

PRINT 1 1 50 ♦ ( X 1 4B ( I ) ♦ I = 1 « 9 ) 
PRINT 1650 

PRINT 1 1 50 » < X I 4B (I)«IslO«l8) 
PRINT 1660 

PRINT 1 150* <XI5B< I)*I=1*9) 
PRINT 1670 

PRINT 1 150* <XI5B< I ) « 1=10* 18) 
PRINT 1680 

PRINT 1 150* (XI6B< I ) » 1=1 *9) 
PRINT 1 690 

PRINT 1150*<XI6B< I )« 1=10*18) 

RETURN 

END 


SLO « • GU 1 * • 0 
$XQ« • MA 1 N ♦ 0 
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SUBROUTINE counts 
data ILOC/O/ 

LOGICAL INDIC 
EXTENDED title 

lM253.nBM"l?^r 0B ’ M25B ’ M90 ''’ M2S ' , ’ M90 ’ M2S,MA25,MA90 ’ MB90 ’ M903 ’ 

^ COMMON BLOCK /- I0CT8. I . IERR.NM4X.NRMAX.NPASS.NRON . INDIC. I CHAN. 

COMMON ADSCAL ( 16) «DASCAL< 16) »BX<4> »DEL4X<3) «DEL5X<3) .DEL fivf 
1 ISX<4> , IM90X.6) , I M 25X<6> , 104(3. . ,05(3. . 1 06 V AB ulw 
2YDDB( 18 ) «PHI2B ( 1 8 > » VW B ( 18 . *ATA1 B( 18 ) « ATA2B (18) t ATA3B ( 18. i ATA4B ( J 8 ) 

= A 3?iS' 1 ® , ’ AT0ZB< 1B).AT03B( 18) » ATD4B, 18) *PH1B( 16. • X I IB ( 18).XI2B( IS) 

;2!“' 1 181 ' 0EU,a ‘ 18 ’ • 0EL2B< 

9AT04, ,8. .«( ,8, .OEL4, ,B, ,DEL5( 18, .DELSUB, .UBX^V. ! ' 

1 I M25X( 6) ♦ I ID4 (3 ) * I 1 05 1 3 ) * I ID6(3 ) .B3( 18 ) . A3< 18) ♦ VW3 ( 18 ) «TI TLE( 20 ) 

o ?^ (18),S33( 18MS34( 18,9541 (18, ‘ S ^< 16). S43 (18). S44 ( 18). 
o , rfp:!« ,,S53<l8>,S54(18)9S6l(18)9S62(18 >^63( 8) ^64 8 

^3^ : is: js; :;f 1 1 is; j; • 

common 

1 M90X ( 6 ) . M25X ( 6 ) . M 90B( 18) .M25B( 18) .M90V< 18) *M25V (18). M90 ( 1 fl 5 • 

1M25< 18) .MA25( 18) .MA90C18) .MB90C 18) .M903( 18) .M253Vl8) 

2116(18). I2B(18).I3BC18) 2 3<18> ’ 

1700 FORMAT ( I 1 4 »8E 1 4. 7 ) 


26 1 I BX ( 1 ) =0 
I IBX( 3) =0 
DO 27 I = 1 ♦ 5.2 
I IM90X( I ) =0 
I I M25X ( I ) =0 
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27 CONTINUE 

C I F < I LOC. AND. 1 ) I IBXI 1 >*1 

C IF ( I LOC* AND *2 ) I 1 BX ( 3 ) s 1 

C I F ( I LOC • AND *4)11 M90X < J ) = 1 

C IF ( I LOC. AND. 8 ) I I M90X( 3 ) * 1 

C I F < I LOC « AND .16) I I M90X ( 5 ) = 1 

C IF< ILOC. AND. 32>I IM25X(1)=1 

C IF( I LOC. AND .64 ) I IM25X<3>*1 

C IF< I LOC. AND. 128) I IM25X<5)=1 

ASSEMBL 

C A /ILOC 

AHM 1 / I IBX+O * • 1 5 

ROT 1 

AHM 1 /I IBX+2. . 15 

ROT 1 

AHM 1 /I IM90X+0* . is 

ROT 1 

AHM1 /I IM90X+2* » 15 

ROT 1 

AHM 1 /I IM90X+4* ♦ 1?$ 

ROT 1 

AHM 1 /I IM25X+0, * 1 «=; 

ROT 1 

AHM 1 /I IM25X+2. » 15 

ROT 1 

AHM 1 /I IM25X+4 * ♦ 1 5 

FORTRAN 

I BX ( 1 ) = I BX ( 1 ) + 1 1 8X ( 1 ) 

I BX ( 3 ) * I BX < 3 ) + 1 I BX ( 3 ) 

DO 28 I*l«5«2 

IM90XI I ) * I M90X ( I )+ I IM90X< I ) 
IM25X ( I ) * I M25X < I >+ 1 1 M25X ( I ) 

28 CONTINUE 

DO 210 1*1.3 
1 ID4( I ) *0 
1 ID5 < 1>*0 
I ID6 ( I ) ~0 
210 CONTINUE 
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DO 200 I = 1 « NMAX 

DEL4 ( I >=.6042*CS41 ( I >*ATA1 { I )+S42( I )*ATA2( I )+S43< I )*ATA3< I ) + 

1 S44( I >*ATA4< I ) ) 

DEL5< I )=.59*(S5l ( I >*A tA1 ( I >+S52< I )*ATA2< I 1+S53I I )*ATA3< I ) 

1 4- S54 ( 1 )*ATA4< I ) ) 

DEL6 ( I ) = «64* ( S61 < I )*ATA1 ( I )+S62( I )*ATA2< I )+S63( I )*ATA3( I )+ 

1 S64( I >*ATA4( I ) ) 

IF ( DEL4 ( I ) *GT *DEL4X ( 1 > ) I ID4( 1 ) = 1 
IF(DEL4( I ) • GT • DEL4X ( 2) ) I ID4(2)*I 
I F ( DEL4 ( I ) .GT.DEL4X(3) ) I I D4 < 3 ) * 1 
I F ( DEL5 ( I ) .GT.DEl5X( 1 ) ) I I D5 < 1 ) * 1 
IF ( DEL5 ( I ) .GT.DEL5X(2) ) I ID5(2)*1 
IF ( DEL5( I ) • GT • DEL5X ( 3 ) ) I ID5( 3) * 1 
IF ( DEL6( I ) *GT*DEL6X( 1 ) ) I ID6( 1 )*I 
I F ( DEL6( I ) »GT *DEL6X C 2 ) > I ID6(2)*1 
I F ( DEL6 ( I ) • GT • DEL6X (3) ) I !D6(3)=1 
200 CONTINUE 

DO 220 1=1*3 

ID4< I )=ID4< I )+I ID4( I ) 

I D5 ( I ) = I D5 ( I ) + 1 I D5 (I ) 

ID6< I )*ID6( I >+ I ID6< I > 

220 CONTINUE 

C CALCULATE THE SUMMATION OF XI 
DO 29 1 = 1* NMAX 
BB< I ) =BB { I )+ B ( I ) 

AB< I >=AB< I ) + A ( I ) 

M90B ( I ) =M90B ( I ) +M90 ( I ) 

M25B ( I ) =M25B { I ) +M25 < I ) 

VWB ( I ) = VWB ( I ) +VW(I) 

ATA1BC I )=ATA1B< I ) +ATa1<I> 

ATA2B < I >=ATA2B< I )+ATA2( I ) 

ATA3B( I )=ATA3B( I )+ ATa3<I> 

AT A4B ( I ) =AT A4B ( I )+ ATo4<I) 

ATD1B( I )=ATD1BU H-ATDl ( I ) 

ATD2B( I ) = ATD2B ( I )+ATD2( I ) 

ATD3B ( I )=ATD3B( I )+ATD3( I > 

ATD4BC 1 > = ATD4B < I >+ATD4< I > 

PH I B < I ) = PH I B ( I ) + PH I ( I 1 


C-23 



LMSC/HREC DI494I0 


YDDB ( I ) » YDDB ( I >+YDD< I ) 

PHI 2B< 1 > =PH 1 2B I I >+PHl2< 1 > 
29 CONTINUE 

return 

END 


C-24 



LMSC/HREC D149410 


SUBROUTINE VAR IAN 
LOGICAL INDIC 
EXTENDED TITLE 

REAL M90X * M25X « M90B » MgSB . M90 V « M25 V « M90 » M25 . MA25 * MA90 * MB90 » M903 * 

1 M2S3 « I 1 B . I 2B » 1 3B 

COMMON / BLOCK / I OCT©* I * I ERR • NMAX • NRM AX * NP ASS * NRUN * INDIC. ICHAn* 

1 T I * T2 *N * ILOC 

COMMON ADSCALt 16) *DASCAL( 16) *BX<4) .DEL4X( 3 ) *DEL5X ( 3 ) «DEL6X(3> » 

1 IBX( 4 ) « I M90X ( 6 > ♦ IM25X (6) *ID4(3)*ID5(3>»ID6(3).BB(18)»AB(18)* 

2YDDB ( 18) .PHI2B( 18) *VWB< 18).ATA1B(18) • ATA2B ( 18>.ATA3B( 18)«ATA4B{ 18) 
3« 

I®? * ATD3B( 18) *■ A TD4B( 18) *PHIB< 18) .XI 1B( 18) *XI2B( 18) 

5 ><I3B( 18) «XI4B< 18) «XI5B( 18) ,XI6B( 18)* DEL 1 B < 18) *DEL2B( 18 ) .DEL3B ( 18 ) 
6.DEL4BC 18) .DEL5B ( 18 > ♦ DEL6B ( 18) * Bv( lQU 

7AV( 18) «VWV< 18) »DAC ( 8) *ADC< 16)*B(18)*A( 18)*YDD( 18) .PHI2< 18) .PHI ( 18) 

fl^f! < i? ),ATA2(I8> ’ ATA3<18) ’ ATA4(18),ATDl( ie > ’ATD2( 18) ♦ ATD3( 18). 
9ATD4( 18 ) . VW( 18) »DEL4( 18 ) *DEL5( 18) «DEL 6 ( 18) • I IBX( 4 ). I IM90X<6) * 

II ?o 1 * I ID4 (3> * 1 ID5<3> * IID6<3 >,B3< 1S).A3( 18).VW3( 18) ♦ TITLE! 20). 

ftIJ/ <18) * Sl2<18> * Sl3<l8),S14<18),S2l(18) * S22 ( 18) *S23( 18) *S24( 18) . 
3S31 (18) *S32< 18) »S33< 18) *S34( 18) *S41 ( 18) »S42( 18)*S43( 18) »S44( 18). 

i?^!!f!’? 5 f n8,,S53<l8, * S54(18) ’ S61(18) * S62<l8) * S 63(18).S64<18), 
®I*J ! 2 ’Iff! 18WT13( ^8) *T14( 18) ♦ T2 1(18)* T22 ( 1 8 ) « T23 (18) * T24( 18), 

8 f! Ilf ’Iff ( 18> * T33< 1 9 ),T34( 18 ) 'P11< 18) *P12(18).P13( 18) *P14( 18). 

i fpiI?!o? 2 f * 185 * P23< 18) ,P24( 18) «P3l ( 18) *P32( 18) *P33( 18) *P34( 18) 

1 * A2 1 B ( 1 8 ) * A22B (18)* A23B (18) 

common 

* M90X ( 6 ) . M25X ( 6 ) * M90B ( 18) .M25B( 18) ,M90V( 18) .M25VC 18) .M90< 18) * 

1M25 ( 18) « MA25 ( 18) *MA90 ( 18) *MB90( 18) *M903( 18) «M253( 18) « 

2I1B(18)»1 2B ( 1 8 ) * 1 3B ( 1 8 ) 

CALCULATE THE VARIANCES 
35 DO 42 1=1. NMAX 

BV ( I ) =8V ( I ) / ( NRMAX— 1 ) 

AV( I )=AV( I ) /'(NRMAX— 1 ) 

M90V ( I )=M90V( I )/ (NRMAx- 1 ) 

M25V ( I ) =M25V ( I ) X ( NRMAx- 1 ) 

VWV ( I ) eVWV ( I )/( NRMAX- 1 ) 

83(1) =3. *SQRT ( BV ( I ) ) 

A3( I ) =3* #SQRT ( AV ( I ) ) 
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M903( I ) =3 • *SQRT ( M90V ( f ) ) 

M253 ( I )*3**SQRT(M25VC I ) ) 

VW3( I )=3.*SQRT< VWV< I ) ) 

XI 1B< I >« SI 1 ( I ) *ATA IB ( I ) +S12< I )*ATA2B< I > +S 13(1) *ATA3B ( I ) 

1 +S 1 4 < I ) *AT A4B ( I ) 

X I2B ( I ) * S21 ( I )*ATA1B( I ) +S22 ( I ) *AT A2B < I ) +S23 ( I ) *ATA3B ( I) 

1 +S24C I)*ATA4B< I ) 

X 1 3B ( I ) x S31 ( I )*ATA1B< I ) +S32 ( I) *ATA2B< I ) + S33 ( I) *ATA3B (I > 

1 +S34( I >*ATA4B{ I) 

X 1 4B ( I ) * S41 < I >*ATAIB< I > +S42 ( I ) *ATA2B( I ) +S43 (I ) *ATA3B ( I > 

I +S44( I )*ATA4B< I) 

X 1 5B < I )= S51 < I )*ATA1B< I > +S52 ( I ) *ATA2B < I ) +S53 ( I ) *ATA3B ( I) 

1 +S54( I >*ATA4B( I > 

XI6BCI)* S61 < I ) *ATA IB ( I ) +S62 ( I ) *ATA2B( I ) +S63 ( I) *ATA3B < I > 

I +S64 ( I )*ATA4B( I ) 

DEL 1 B < I > = XI 1B< I >*,5229 
DEL2B < I )= XI2B(I)**539 
DEL3B( I >= XI3B<I)**543l 
DEL4B( I > = XI4B( I >**6042 
DEL5B(I>= XI5B< I )*»59 
DEL6B < I ) = X 1 6B ( I ) * • 64 

A21B( I )*-YDDB( I )+ T 1 1 ( I >*ATD1B< I >+T12< I )*ATD2B( I ) +T13( I )*ATD3B( I ) 
1 + T 14 ( I )*ATD4BC I ) 

A22B < I )=-YDDB< I )+ T21 ( I >*ATD1B( I >+T22( I )*ATD2B( I ) +T23( I >*ATD3B< I ) 
1 4- T14 ( I ) *ATD4B< I) 

A23B ( I ) =— YDDB ( I )+ T31 ( I )*ATD1B( I )+T32< I )*ATD2B( 1 ) +T33I I )*ATD3B< 1 1 
1 + T34 ( I )*ATD4B< I ) 

I 1B( 1) =~PH I 2B ( I ) +P1 l { I )*ATD1B< I >+Pl2< I )*ATD2B( I ) +P 1 3 ( I >*ATD3B< I ) 
1 +pi4< I )*ATD4B( I ) 

1 2B ( I ) =-PHI28< I > +P2 1 ( I )*ATD1B< I > +P22 ( I >*ATD2B< I >+P23< I )*ATD3B( I > 
1 +P24I I )*ATD4B( I ) 

I 3B ( I ) = -PHI2B( I ) +P31 ( I )*ATD1B( I ) +P32 ( I >*ATD2B< I )+P33< I )*ATD3B( I ) 
1 +P34< I)*ATD4B< I) 

42 CONTINUE 
RETURN 
END 
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subroutine read 

REAL M90X « M25X . M90B • M25B . M90V * M25V • M90 * M25 * MA25 . M A90 . MB90 ♦ M903 . 
1M253. 1 1 B. I 2B. I3B 

COMMON ADSCAL < 16 > »DASCAL( 16) ♦BXt 4 ) *DEL4X ( 3 ) .DEL5XI3) «DEL 6 X( 3) » 

1 I BX < 4 ) ♦ I M90X ( 6 ) * IM25X (6) . ID4(3) . ID5(3> . ID6<3) .BBC 18) .ABC 18 ) ♦ 

2YDDBC 18) .PHI2BC 18) .VWbC 18) . ATA1BC 18) . ATA2BC 18) . ATA3BC 18) . ATA4B C 18) 
3. 

4ATD1BC 18) . ATD2BC 18 ) .ATD3BC 18) .ATD4BC 18) .PHIBC 1 8 ) . X 11 B C 1 8 ) . X I 2B C 18 ) 

5. XI3BC 18) .XI4BC 18) .XIsBC 18) .XI 6 BC 18) . DEL 1 B ( 18) .DEL2BC 18 ) .DEL3BC 18) 

6. DEL4BC 18) .DEL5B C 18) »DEL6BC 18) « BV<18>* 

7 AVC 18) .VWVC 18) .DACC8) «ADC( 16)«BC 18). AC 18) . YDD ( 18) .PHI2C 18) .PHI C 18) 
8 » AT A 1 ( 18) ♦ ATA2C 18) . ATA3C 18) .ATA4C 18) « ATD 1 C 18) ♦ ATD2C 18) * ATD3C 18) . 
9ATD4C 18) .VWC 18) .DEL 4 C 18 ) . DEL 5 C 18) « DEL 6 C 18) . I IBXC4 ) . I IM90XC6) ♦ 

1 I IM25XC 6 ) . I ID4C3 > ♦ I 105(3) « I ID6C3) »B3< 18) .A3( 18) *VW3( 18) *TITLE(20) . 
2S11 C 18) .S12C 18).S13C 18) «S14< 18) .S21 C 18) .S22C 18).S23C 18) .S24C 18). 
3S31 C 18) .S32C 18) .S33C 18) .S34C 18) .S41 ( 18) .S42( 18).S43< 18) .S44C 18). 
4S51 C 18) » S52 C 18) .S53C 18) «S54C 18) .S61 C 18) *S62( 18) »S63C 18) .S64C 18 ) ♦ 

5T 11(18). T12(18).T13(18)*T14(18).T2 1(18) .T22( 18) »T23( 18) *T24( 18) . 
6T31 ( 18) .T32C 18 ) .T33( 18) .T34( 18).Pll(18).P12(18).P13(18).Pl4(l8). 
7P21 ( 18 ) .P 22 ( 18 ) .P23C 18) »P24( 18) «P31 ( 18 ) .P32( 18 ) .P 33 ( 18 ) .P 34 ( 18 ) 

1 . A 21 BC IT ) . A 22 BC 18 ) ♦ A 23 B ( 18 ) 

COMMON 

1 M90X(6) . M25X ( 6 ) .M90b( 18) .M25B( 18) .M90V( 18) .M25V( 18) .M90( 18) . 

1M25 ( 18) » MA25 ( 18 ). MA90 ( 18 ). MB90 ( 18 ). M903 (18). M253 ( 18 ). 

21 1B( 18) . I2B( 18) * I3B( 18) 

1020 FORMATC 8 EIO. 4 ) 

C INPUT A/d AND D/A SCALE FACTORS 
READ (5*1 020 ) ADSCAL 
READ (5.1 020 ) DASCAL 
C INPUT EXCEEDANDE LEVELS 
READ (5. 1020) BX 
READ (5* 1020) M90X 
READ C 5 . 1020) M25X 
READ (5. 1020) DEL4X 
READ ( 5 * 1020) DEL5X 
READ (5* 1020) DEL 6 X 
READ ( 5 . 1020) MA25 
READ (5* 1020 ) MA90 
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READ < 5 « 1 020 ) MB90 
READ <5. 1020) Sll 
READ (5* 1020) s 12 
READ ( 5* 1 020 ) S13 
READ <5* 1020) S14 
READ <5* 1020) S21 
READ (5* 1020) S22 
READ (5*1020) S23 
READ (5* 1020) S24 
READ <5* 1020) S31 
READ <5* 1020) S32 
READ (5* 1020) S33 
READ (5* 1020) S34 
READ <5*1 020 ) S41 
READ (5* 1020) S42 
READ (5* 1020) S43 
READ (5* 1020) S44 
READ (5* 1020) S51 
READ (5* 1020) S52 
READ (5* 1020) S53 
READ (5* 1020) S54 
READ (5* 1020) S61 
READ(5«1020) S62 
READ (5* 1020) S63 
READ (5* 1020) S64 
READ (5* 1020) Til 
READ (5*1 020 ) TI2 
READ <5* 1020) T13 
READ (5* 1020) T 14 
READ (5* 1020) T21 
READ <5* 1020) T22 
READ (5* 1020) T23 
READ <5* 1020) T24 
READ <5* 1020 ) T31 
READ (5* 1020) T32 
READ (5* 1020) T33 
READ (5* 1020) T34 
READ (5*1 020 ) Pll 
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READ ( 5» 1020) 
READ< 5* 1020) 
READ (5* 1020) 
READ( 5« 1020) 
READ ( 5« 1020) 
READ! 5* 1020) 
READ <5* 1020) 
READ! 5* 1020) 
READ (5* 1020) 
READ* 5* 1020) 
READ (5* 1020) 

return 

END 


P12 
P13 
PI 4 
P21 
P22 
P23 
P24 
P31 
P32 
P33 
P34 
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